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Abstract 

In the leading twist approximation of the Wilson operator product expansion with 
standard and "frozen" versions of strong copling constant we show that the Bessel- 
inspired behavior of the structure function F% at small x, obtained for a flat initial 
condition in the DGLAP evolution, leads to very good agreement with the deep 
inelastic scattering experimental data from HERA. 



1 Introduction 

The measurements of the deep-inelastic scattering structure function (SF) F 2 at HERA pQ 
have permitted the access to a very interesting kinematical range for testing the theoretical 
ideas on the behavior of quarks and gluons carrying a very low fraction of proton momentum, 
the so-called small x region. In this limit one expects that the non-perturbative effects may 
give an essential contributions. However, the reasonable agreement between the HERA data 
and the next-to-leading order (NLO) approximation of perturbative QCD, which has been 
observed for Q 2 > lGeV 2 (see the reviews in [2]), indicates that the perturbative QCD could 
describe the SF evolution up to very low Q 2 values, traditionally explained by soft processes. 
It is of fundamental importance to find out the kinematical region where the well-established 
perturbative QCD formalism can be safely applied at small x. 

The standard program to study the small x behavior of quarks and gluons is carried out 
by comparison of the data with the numerical solution of the DGLAP equations fitting the 
parameters of the x profile of partons at some initial Ql and the QCD energy scale A (see, 
for instance, [31 H]). However, in analyzing exclusively the small x region (x < 0.01), there is 
the alternative of doing a simpler analysis by using some of the existing analytical solutions 
of DGLAP in the small x limit (see [2] for review). It was done in Refs. [5]- [8], where it 
was pointed out that the HERA small x data can be interpreted in the so called doubled 
asymptotic scaling (DAS) approximation related to the asymptotic behavior of the DGLAP 
evolution discovered in [9] many years ago. 

Here we illustrate results obtained in [6j-[8j, which are the extension of previous leading 
order (LO) studies [9l [5] to the NLO QCD approximation. The main ingredients are: 

1. Both, the gluon and quark singlet densities are presented in terms of two compo- 
nents ('+' and '— ') which are obtained from the analytical Q 2 dependent expressions of the 
corresponding ('+' and '— ') parton distribution (PD) moments. 

2. The '— ' component is constant at small x, whereas the '+' component grows at Q 2 > Qq 
as ~ exp (cryvLo)) where 



°nlo = 2\J (d + s + D + p) \n(x), 



the LO term d+ = — 12//3 and the NLO one D+ = d ++ + d+fa/Po wi th d ++ = 412//(27/3b). 
Here the coupling constant a s = a a /(47r), s = Zn[a s (<5o)/a s (<3 2 )] an d p = a s (<3o) — a s(Q 2 )> Po 
and /3i are the first two coefficients of QCD /^-function and / is the number of active flavors. 



2 Basic formulae 

Our purpose is to extract the small x asymptotic PD form in the framework of the DGLAP 
equation starting at some Qq with the flat function: 

f T a \Ql) = A a (hereafter a = q,g), (1) 

where f r2 are the leading-twist parts of parton (quark and gluon) distributions, multiplied 
by x, and A a are unknown parameters that have to be determined from data. We neglect the 
non-singlet quark component at small x. 

We would like to note that HERA data [lj show a rise of F 2 at low Q 2 values (Q 2 < lGeV 2 ) 
when x — > (see Fig.l below). This rise can be explained naturally by incorporation of higher- 
twist terms in the analysis (see [8] and Fig.l). 

We shortly compile below the LO results (the NLO results may be found in [6J [8]). The 
full small x asymptotic results for PD and SF F 2 at LO are: 

F 2 (x,Q 2 ) = ef q (x,Q 2 ), (2) 
f a (x,Q 2 ) = f+(x,Q 2 ) + f-(x,Q 2 ) , (3) 

where e = (Xa e 2 )// is the average charge square. The '+' and '— ' components fai x , Q 2 ) are 
given by the sum 

f?(*, Q 2 ) = f?*(x, Q 2 ) + fa T 'H^ Q 2 ) (4) 

of the leading-twist parts f^ 2,± (x,Q 2 ) and the higher-twist parts f^ T,± (x,Q 2 ), respectively. 
The small x asymptotic results for PD f^ 2,± are 

/f' + M 2 ) = (A g + ^A q y o (a)e-~ d +Ws + 0(p) , (5) 

r 2 > + (x, Q 2 ) = f -(A g + U q )p h{a) e-~ d ^> + 0{p) , (6) 
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f^-(x,Q 2 ) = ~A q e- d -^ s + 0(x), (7) 

f^-(x,Q 2 ) = A q e' d -^ s + 0(x), (8) 

where d+(l) = 1 + 20//(27/3 ) and cf_(l) = 16//(27/3 ) are the regular parts of d+ and oL 
anomalous dimensions, respectively, in the limit n — > 1 0. The functions l v {y = 0, 1) are 
related to the modified Bessel function I v and to the Bessel function J u by: 

h[a) ~ \ J u (a), if s < ■ (9) 
The variables a and p are given by 



a = 2Vrf +S ln(x), P=^ = ^^y (10) 



1 For a quantity k{n) we use the notation k{n) for the singular part when n — + 1 and k(n) for the 
corresponding regular part. 



3 Effective slopes 



As it has been shown in Refs.[6]-[8], the PD and F 2 behavior, given in the Bessel-like form 
by generalized DAS approach, can mimic a power law shape over a limited region of x and 
Q 2 : 

f a (x, Q 2 ) ~ x - xl "^ Q2) and F 2 (x,Q 2 )~ x - x r"^ 2 \ 
At the twist-two LO approximation, they have the following form 

fg{X,Q 2 ) Iq(<t) 

\eff (r q2\ _ X eff( Q 2\ _ fg( X iQ ) hi?) (n) 

X F2 (x,Q) - X q [x,Q) - fq{X)Q 2 ) P ~ h{a y I 11 ) 

The corresponding NLO expressions and the higher-twist terms can be found in Refs. [6]- [8]. 

The effective slopes \ e Jf and X e J 2 depend on the magnitudes A a of the initial PD and 
also on the chosen input values of Ql and A. To compare with the experimental data it 
is necessary to use the exact expressions ( 1TTT) but for qualitative analysis one can use some 
appropriate approximations. 

At large values of Q 2 , the "— " component of PD is negligable and the dependence of slopes 
on the PD disappears. In this case the asymptotic behaviors of slopes have the following 
expressions 0: 

K ff ' as &Q 2 ) = ptQ*p 



f (o) 41n(l/x)' 

x,Q 2 ) = Xf' as (x,Q 2 ) = P £M„ P _ 3_ ( 12 ) 

I Jo) 41n(l/x) 



where the symbol ~ marks the approximation obtained by the expansion of the usual and 
modified Bessel functions in (jUj). 

One can observe from (|T2|) . that the gluon effective slope \ e Jf ,as is larger than the quark 
slope \gff ,as , which is in excellent agreement with global analyses [3J. 



4 Comparison with experimental data 

With the help of the results presented in the previous sections we have analyzed F 2 HERA 
data at small x from the HI and ZEUS collaborations pQ. In order to keep the analysis as 
simple as possible we have fixed the number of active flavors f=4 and ^j^g( n / = 4) = 292 
MeV in agreement with the more recent HI results |10j . 

The typical fits for the SF F 2 (x, Q 2 ) as a function of x for different Q 2 bins are presented 
in Fig. 1. The experimental points are from HI (open points) and ZEUS (solid points). 
The solid line represents the NLO fit alone with x 2 /n.d.f. = 1.31. The dashed curve are 
obtained from the fit at the NLO, when the renormalon contributions of higher-twist terms 
have been incorporated. The corresponding x 2 /n.d.f. = 0.86. The dash-dotted curve (hardly 

2 The asymptotic formulae given in Eq. (|12p work quite well at any Q 2 > Qg values, because at Q 2 = Q\ 
the values of \ e J^ and are equal zero. The use of approximations in Eq. (p~2|) instead of the exact results 
given in Eq. (flTj) underestimates (overestimates) only slightly the gluon (quark) slope at Q 2 > Qq. 



distinguished from the dashed one) represents the fit at the LO together with the renor- 
malon contributions of higher-twist terms. The corresponding x 2 / n -d-f- = 0.84. The results 
demonstrate excellent agreement between theoretical predictions and experimental data for 
the region Q 2 > 0.5 GeV 2 . However, the twist-two appriximation is in agreement with the 
data only for Q 2 > 2.5 GeV 2 . 

Using these results of the fits of the SF F 2 (x, Q 2 ) we analyze also the HERA data for the 
slope din i^/cf ln(l/x) at small x from the HI and ZEUS Collaborations [T0]-[T2]. The results 
are shown in Fig. 2. Because the twist-two approximation is reasonable at Q 2 > 2.5 GeV 2 
only, some modification should be considered at the lower Q 2 values. In the paper [8] we have 
added the higher twist corrections and have found a good agreement for Q 2 > 0.5 GeV 2 . 

Here we study another possibility. We modify the QCD coupling constant. From different 
studies [T3"t [14] it is known that the effective argument of the coupling constant is higher then 
Q 2 at low x values. In the present paper we consider freezing of the coupling constant by 
changing its argument Q 2 —> Q 2 + M 2 , where M p is the p- meson mass [T5]. Thus, in the 
formulae of the Section 2, we should do the following replacement 

a s (Q 2 ) -> a fr (Q 2 ) = a s (Q 2 + M 2 ) (13) 

Fig. 2 shows the experimental data for \p 2 . It demonstrates an essential improvement 
between the low Q 2 ZEUS data for A"^//(x, Q 2 ) and our results with the "frozen" coupling 
constant df r {Q 2 ). 

Indeed, from our fits for F 2 (x, Q 2 ) we have found that Ql ~ 0.5-T-0.8 GeV 2 j8]. So, initially 
we had X e p 2 (x, Ql) = from the suggestion ([I]). The replacement ({TBI modifies the value of 
X e p2 {x,Ql). For the "frozen" coupling constant aj r (Q 2 ) the value of \p 2 {x,Ql) is non-zero 
and the slope is quite close to experimental data at Q 2 ~ 0.5 GeV 2 . Nevertheless, for Q 2 < 0.5 
GeV 2 there is still an disagreement with the data, that needs additional investigations. 

5 Conclusions 

We have shown that the results, developed recently in [5]- [8], have quite simple form and 
reproduce many PD properties at small x, that have been known from global fits. 

We found the very good agreement between our approach, based on QCD, and HERA data, 
as it has been observed earlier with other approaches (see reviews [2]). The (renormalon-type) 
higher-twist terms lead to the natural explanation of the rise of SF F 2 at low values of Q 2 
and x, which has been discovered in recent HERA experiments [JJ. 

The application of the "frozen" coupling constant a,f r (Q 2 ) leads to good agreement with 
the recent HERA data PJj]-[I2] for the slope X e J 2 f (x,Q 2 ) for Q 2 > 0.5 GeV 2 . As the next 
step of our investigations, we plan to fit the HERA experimental data for the F 2 (x,Q 2 ) SF 
directly with the "frozen" coupling constant af r (Q 2 ). 
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